Transition from localized to extended eigenstates in the ensemble 



of power-law random banded matrices 

Alexander D. Mirlin^f, Yan V. Fyodorov^f, 
Frank-Michael Dittes^, Javier Quezada^ and Thomas H. Seligman^ 

■ ^ Institut fiir Theorie der Kondensierten Materie, Universitdt Karlsruhe, 76128 Karlsruhe, 



p) • Fachbereich Physik, Universitdt-GH Essen, 45117 Essen, Germany 

^ Forschungszentrum Rossendorf, Institut fiir Kern- und Hadronenphysik, 01314 Dresden, 



(N 



Germany 



Germany 



^ . ^ Tecnologico de Monterey, Guadalajara Gampus, Guadalajara, Mexico 

. ^ IFUNAM Laboratorio de Guernavaca, 62191 Guernavaca, Morelos, Mexico 

5 . (February 1, 2008) 

S Abstract 

\ We study statistical properties of the ensemble of large N x N random 

Q I matrices whose entries Hij decrease in a power-law fashion Hij ~ \i — 

J-^ ' Mapping the problem onto a nonlinear cr— model with non-local interaction, 

we find a transition from localized to extended states at a = 1. At this 

^ . critical value of a the system exhibits multifractality and spectral statistics 

^ i intermediate between the Wigner-Dyson and Poisson one. These features are 

reminiscent of those typical for the mobility edge of disordered conductors. We 

find a continuous set of critical theories at q = 1, parametrized by the value 

of the coupling constant of the cr— model. At a > 1 all states are expected to 

be localized with integrable power-law tails. At the same time, for 1 < a < 

3/2 the wave packet spreading at short time scale is superdiffusive: (|r|) ~ 
1 

t2Q-i, which leads to a modification of the Altshuler-Shklovskii behavior of 
the spectral correlation function. Atl/2<Q!<1 the statistical properties 
of eigenstates are similar to those in a metallic sample in d = {a — l/2)~^ 
dimensions. Finally, the region a < 1/2 is equivalent to the corresponding 
Gaussian ensemble of random matrices (a = 0). The theoretical predictions 
are compared with results of numerical simulations. 
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I. INTRODUCTION. 



Recently, there has been a considerable interest in the properties of large N x N ran- 
dom banded matrices (RBM). The ensemble of RBM is defined as the set of matrices with 
elements 

H,, = G,,a{\t - j\) , (1) 

where the matrix G runs over the Gaussian Orthogonal Ensemble (GOE), and a(r) is some 
function satisfying the condition limj.^00 cii'f') = and determining the shape of the band. 
In the most frequently considered case of RBM the function a(r) is considered to be fast (at 
least, exponentially) decaying when r exceeds some typical value b called the bandwidth. 
Matrices of this sort were first introduced as an attempt to describe an intermediate level 
statistics for Hamiltonian systems in a transitional regime between complete integrability 
and fully developed chaos [|l| and then appeared in various contexts ranging from atomic 
physics (see P] and references therein) to solid state physics P] and especially in the course 
of investigations of the quantum behavior of periodically driven Hamiltonian systems [Q. 
The mostly studied system of the latter type is the so-called quantum kicked rotator (KR) 
^ characterized by the Hamiltonian 

j2 00 

H = ^ + V{9) Y: 5it-mT), (2) 

m=—oo 

where / = —ihd/dO is the angular momentum operator conjugated to the angle 9. The 
constants T and I are the period of kicks and the moment of inertia, correspondingly, and 
V{9) is usually taken to be V{9) = k cos 9. Classically, the KR exhibits an unbound diffusion 
in the angular momentum space when the strength of kicks k exceeds some critical value. It 
was observed, however, that in a quasi-classical regime quantum effects suppress the classical 
diffusion [^] in close analogy with the effect of Anderson localization of a quantum particle 
by a random potential [^]. 

It is natural to consider the evolution (Floquet) operator U that relates values of the 
wavefunction over one period of perturbation, ipi^^t + T) = Uip{9,t), in the "unperturbed" 
basis of eigenfunctions of the operator /: |/) = (^2^)^/2 exp(m6'), n = ±0, ±1, . . .. The matrix 
elements {m\U\n) tend to zero when |m — n| — 00. In the case V{9) = kcos9 this decay 
is faster than exponential when |m — n\ exceeds b ~ k/h, whereas within the band of the 
effective width b matrix elements prove to be pseudorandom ||^. All these observations 
attracted much interest to the statistical properties of the ensemble of RBM in order to use 
the extracted information for understanding the dynamical properties of the KR. 

Let us note, however, that the fast decay of (m|[/|n) in the above mentioned situation 
is due to the infinite differentiability of V{9) = cos9. If we took a function V{9) having a 
discontinuity in a derivative of some order the corresponding matrix elements of the evo- 
lution operator would decay in a power-law fashion when \n — m\ oo[|. In fact, there 



The authors are grateful to F. Izrailev for attracting their attention to this fact. 
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is an interesting example of a periodically driven system where the matrix elements of the 
evolution operator decay in a power-law way, namely the so-called Fermi accelerator 0. 
This system does not show the effect of dynamical localization in energy space typical for 
the KR. One can hope to understand this difference in behavior studying properties of RBM 
with power-law decay of the function a(r) at infinity. 

One may also consider the random matrix (1) as the Hamiltonian of a one-dimensional 
tight-binding model with long-ranged off-diagonal disorder (random hopping). A closely 
related problem with non-random long-range hopping and diagonal disorder was studied 
numerically in Ref. [Q. The qualitative effect of weak long-range hopping on the localized 
states in a 3D Anderson insulator was discussed by Levitov 0. We will discuss the corre- 
spondence between our results and those of Ref. later on. Let us note that similar models 
with a power-law hopping appear also in other physical contexts ||1CI|| . 



As it was shown in Ref. the conventional RBM model can be mapped onto a ID 

supermatrix non-linear a— model, which allows for an exact analytical solution. The same 
cr— model was derived initially for a particle moving in a quasi- ID system (a wire) and being 
subject to a random potential. All states are found to be asymptotically localized, with the 
localization length 



^ = ^-^2 0^ ^ ' Bk= a (r)r . (3) 

^-DQ r=-oo 



In the present paper we consider the case of a power-like shape of the band: 

a(r) oc r~" for large r . (4) 



Under these conditions, the derivation of the cr— model presented in [^,rT| loses its validity. 



Reconsidering this derivation, we arrive at a more general ID a— model with long-range 
interaction. Performing its perturbative and renormalization group analysis, we find that 
the model is much richer than the conventional short-range one. In particular, it exhibits an 
Anderson localization transition at a = 1. The main scope of the present work is to study 
the statistical properties of the model in the whole range of the parameter a. 



II. THE POWER-LAW RANDOM BANDED MATRIX ENSEMBLE AND THE 

EFFECTIVE NON-LINEAR cr-MODEL 

Let us consider the ensemble given by eq. (|I|) with the function a(r) having the form 

«W = I (r/fe)-", ~r>h ■ 

The parameter h will serve to label the critical models with a = 1. We will consider h to 
be large: 6 ^ 1, in order to justify formally the derivation of the cr— model. We will argue 
later on that our conclusions are qualitatively valid for arbitrary h as well. We will call the 
ensemble ([l|), (|^) the power-law random banded matrix (PRBM) model. 

In Ref. [0,0 it was shown that the RBM model (|ip can be mapped, for arbitrary 
bandshape a(r), onto a field-theoretical model of interacting 8x8 supermatrices cr^ , = 
1,2,..., A^) characterized by the action 
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S{(j} = Str J i ^ RijCTiCTj + YlUicri 
Here Rij = (A''^)ij — -^6ij Y.ki{^^^)kh ^ is the matrix with elements Aij = a^{\i — 



(6) 



U{a) = istr hn{E~a- ^^A) + 1 Ili^-'h^' 



(7) 



CO is the frequency, and A = diag(J, —I) in the "advanced-retarded" representation. The 
supermatrix cxj can be parametrized as cTj = Tj^^PiTi, with Pi being block- diagonal in the 
"advanced-retarded" decomposition and Tj belonging to a certain graded coset space, see 
the reviews |T^J13[| . For 6^1 the integral over the matrices Pi can be evaluated by the 
saddle-point method |TT[. Then the action is reduced to a cr— model on a lattice: 



(8) 



Here Qi = Ti ^ATj satisfies the constraint = 1, Aq is given hyf\ Aq = J2i Aki 
I]J^_oo z/ is the density of states: 



2nA 



1/2 



(9) 



The standard next step is to restrict oneself to the long wavelength fiuctuations of the 
Q— field. For usual REM characterized by a function a{r) decreasing faster than any power 
of r as r — i> oo, this is achieved by the momentum expansion of the first term in the action 



1] ^0 



QiQj = E [^9 ^ ~ ^0 ^ QqQ 



1-3 



B2_ 

2Ao 



2An 



dx{d,Qf 



(10) 



where B2 = J2i^ki{k — /)^, as defined in eq. (^. This immediately leads to the standard 
continuous version of the nonlinear cr— model: 



S{Q} 



Str / dx 



^Do {d^Qf + iujQK 



(11) 



with the classical diffusion constant Dq = 7rz/i?2, which implies the exponential localization 
of eigenstates with the localization length C, = ttuDq oc 6^. 



^The expression for Aq is valid for a > 1/2 and in the limit N — > 00. When a < 1/2, Aq starts to 
depend on N which can be removed by a proper rescaling of the matrix elements Hij in eq. (|l|). 
Then, the properties of the model with a < 1/2 turn out to be equivalent to those of the GOE, so 
we will not consider this case any longer. 
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Let us try to implement the same procedure for the present case of power-hke bandshape, 
eq. (^. Restricting ourselves to the lowest order term in the momentum expansion, one 
arrives again at eq. ( |Tl| ) as long as a > 3/2. This suggests that for a > 3/2 the eigenstates 
of the present model should be localized in the spatial domain of the extension ^ oc i'Dq. 
However, in contrast to the usual REM model we expect this localization to be power-like 
rather than exponential: |'?/'(r)p oc r~^" at r ^ This is quite evident due to the possibility 
of direct hopping with the same power-law. On a more formal level the appearance of power- 
law tails of wavefunctions is a consequence of the breakdown of the momentum expansion 
for the function — Aq^ in higher orders in g^. The presence of power-law "tails" of 
the wave functions, with an exponent a determined by the decay of hopping elements, was 
found in numerical simulations in Ref. 

The most interesting region l/2<a<3/2 requires a separate consideration. The matter 
is that eq. (|10D loses its validity in view of the divergency of the coefficient B2. Instead, a 
close inspection shows that 



Al{A-' - A^') ^Ao-Aq = 2j^ dra\r){l - cosgr) 

y^^'^ dx{l - cos x) + {b\q\y" j^^ 



2 f r^\i\ , . ,.0., /"^ dx ,1 

— (1 - cosx) ) 

q\ [Jo ' ■ ' Jh\q\ 'J 

Ca&^"|gr°"^ for 1/2 < a < 3/2 and |g| < 1/6 , (12) 



where = 2 "~ cosx) is a numerical constant. 

The corresponding long wavelength part of the action, 

So{Q} = -^Str j rfg|g|2-iQ,Q_, , (13) 

can not be reduced to the local-in-space form in the coordinate representation any longer. 
Here, 1/t = jiiipY Cah'^" oc 6^"~^ ^ 1 plays the role of coupling constant, justifying the 
perturbative/renormalization group treatment of the model described in the next sections. 

Let us mention, that considering the REM model as a tight-binding Hamiltonian, the 
corresponding classical motion described by the master equation on the same ID lattice is 
super diffusive for 1/2 < a < 3/2: (|r|) oc ^^/(^"-i). As will be discussed in Sec. 0, this 
influences the asymptotic behavior of the spectral correlation function for the corresponding 
quantum system. 



III. PERTURBATIVE TREATMENT OF THE NON-LINEAR a-MODEL: 

GENERAL FORMULAS. 

In this section, we derive one-loop perturbative corrections to the density-density corre- 
lation function and inverse participation ratios. The analysis of these expressions for various 
values of the power-law parameter a will be presented in Sec. 0. 
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A. Density density correlation function. 



The basic object characterizing the behavior of a particle in a random medium is the 
density-density correlation function, which can be generally expressed in terms of the 
cr— model as follows ||12|| : 

K{n,r2; uj) = -(7rz/)2 J DQ Qu,a(skf3pQ2i,/3ae-''^'^^ • (14) 

Here the indices p,p' of the matrix Qppi^ap correspond to its advanced-retarded block struc- 
ture, whereas discriminate between bosonic and fermionic degrees of freedom. The 
matrix kjsp is equal to 1 for bosons and (-1) for fermions. In order to calculate the corre- 
lation function (|T^ perturbatively, one needs to parametrize the matrix Q in terms of the 
independent degrees of freedom. We find it convenient to use the following parametrization 
Hi: 

Q = K{W + VlTW^)=x{l + W + — , (15) 

where W is block-off-diagonal in the advanced-retarded representation. To get the pertur- 
bative expansion for K{ri,r2;uj), one has to substitute eq. (p!5|) into (p!^, to separate the 
part quadratic in W from the rest in the exponent and to apply the Wick theorem (see Ref. 
Q for the contraction rules). In the usual case, when the action is given by eq. ([TI|), the 



leading order (tree level) result reads in momentum space as follows: 

Ko{q.^)= „ ■ • (16) 

The perturbative quantum corrections do not modify the general form (p!6D, but change the 
value of the diffusion constant. In particular, in one- loop order one gets eq. (|16D with Dq 
replaced by [|l^ 



D = DAI- —— y } . (17) 

y qi=TTni/Li ) 

This induces the standard weak-localization correction to the conductivity. 

Now we implement an analogous procedure for the non-local a— model of the type of eq. 

(0): 

5{g} = ^Str^f/(r-r')Q(r)Q(r')-2^EStrAg(r) , (18) 

with the Fourier-transform of U{r) behaving at small momenta as 

U[q) = -Igl'^; 1/2 < a < 2 . (19) 

The exponent a is related to the parameter a of the REM model by o" = 2a — 1. In leading 
order, we keep in the action the terms quadratic in W only, which yields: 
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^o(g, uj) = ^"l!" . , (20) 

corresponding to a super diffusive behavior. 

To calculate the one- loop correction to KQ{q) (we set = for simplicity) we expand 
the kinetic term in S{Q} up to fourth order in W: 

Y.^iTU{r-r')Q{r)Q{r') Uth order= E ^Str H^2(r)1^2(r') . (21) 

r,r' r,r' 

The contraction rules are given by eqs. (8), (16) of Ref. with the propagator Il{q) 
replaced by 

n(.) - 4 . (22) 

There is only one one-loop diagram contributing to the self-energy part in the present 
parametrization, see Fig. 1. Evaluating it, we find: 

S^? = Uidk)^Vr^^ (23) 



2J' ' \k\' 

where J{dk) = / 1^ = NJ2k- There is also a contribution to F^^^ from the Jacobian of 
transformation ([15|), which is equal to |1^ 



6T 



(2) 



-^6{0) ^ -IJ (dk) . (24) 



Combining eqs. (E3) and (p4l), we get 



Finally, we get the following expression for the density- density correlation function up 
to one-loop order: 

K~\q) = K^\q) - ^ j {dkY-^ . (26) 



B. Inverse participation ratio. 

In order to characterize eigenfunctions statistics quantitatively, it is convenient to in- 
troduce a set of moments Iq = J2 \'ip{r)\'^'^ of the eigenfunction local intensity |'?/'(r)p. The 
quantity I2 is known as the inverse participation ratio (IPR). The whole set of moments Iq 
is a useful measure of the eigenfunction structure. For completely "ergodic" eigenfunctions 
covering randomly but uniformly the whole sample, (Iq) = as in GOE. In contrast, 

if eigenfunctions are localized in a domain of size ^, then (Iq) oc 1/^"^"^. Finally, the multi- 
fractal structure of the wavefunction manifests itself via the dependence (Iq) oc 1/A^'''j('?-i) 
with dq < 1 being the set of fractal dimensions |jl^,|l8|. 
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The method allowing one to calculate perturbative corrections to the GOE-like results in 
the weak localization regime was developed in 1 1^,14]. It is straightforwardly applicable to 
the present case of power-law RBM, provided the appropriate modification of the diffusion 
propagator entering the contraction rules is made, see the text preceding eq. (|^). One finds 



(27) 



where n(r, r') = (1/A^) J^q^io.) exp[ig(r — r')] and n(g) is given by eq. (^21) . 



IV. RENORMALIZATION GROUP TREATMENT. 



Our effective a— model, eq. ([T8|), is actually of one-dimensional nature. However, for the 
sake of generality, in the present section we find it convenient to consider it to be defined in 
(i-dimensional space with arbitrary d. The form (^) of the generalized diffusion propagator 
suggests that d = a should play the role of the logarithmic dimension for the problem. In the 
vicinity of this critical value it is natural to carry out a renormalization group (RG) treatment 
of the model. We will follow the procedure developed for general non-linear cr— models in 
T^ . We start from expressing the action in terms of the renormalized coupling constant 



t = Zi ^tBl-i''^ where ts is the bare coupling constant and /i ^ is the length scale governing 
the renormalization 0: 



5 = ^^^f/(|r-r |)Str -W{r)W{r') + ^1 + W^ir)^! + W^{r') 



tTCUUJ 



(28) 



Expanding the action in powers of W{r) and keeping terms up to 4-th order, we get 



d—cr 



Ui\r - r'\) Str {W{r) - W{r')f 



AtZi ^ 



d—a 



Si = ^EU{\r- r'DStr W'{r)W'{r') + ^ J] Str W'{r) 

1 7^1-/ T 



(29) 



We have restricted ourselves to 4-th order terms, since they are sufficient for obtaining the 
renormalized quadratic part of the action in one-loop order. The calculation yields, after 
the cancellation of an J{dk) oc 5(0) term with the contribution of the Jacobian: 



,d~<j 



{dk) 



U{k) -U{k + q) 
-U{k)-'-^ 



(30) 



^Note that the VF-field renormalization is absent due to the super symmetric character of the 



problem, which is physically related to the particle number conservation [12 
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According to the renormalization group idea, one has to chose the constant Zi{t) = l+at+. . . 
so as to cancel the divergency in the coefficient in front of the leading Igl*^ term. This will 
be done in the next section. 



V. ANALYSIS OF THE MODEL FOR DIFFERENT VALUES OF THE 

EXPONENT 0- = 2a - 1. 

A. l<cj<2(l<a< 3/2). 

To evaluate the one-loop correction (^) to the diffusion propagator, we use the expansion 



\q + k 


a 




k 


a 



P 2 P V2 MP 



(31) 



. 1^1" ' 

Thus, the integral in eq. (^) can be estimated as 

For a > d = 1 the integral diverges at low k, and the second term in eq. (|32|) dominates. 
This gives 

J~L"-^|g|", (33) 

where L is the system size determining the infrared cut-off (in the original RBM formulation 
it is just the matrix size N). This leads to the following one- loop expression for the diffusion 
correlator 

K{q) = ^^^^ ; r 1 = t-^ - const L''-^ . (34) 

Now we turn to the renormalization group analysis, as described in Sec. |^ In one- 
loop order, the expression for the renormalization constant Zi is determined essentially by 
the same integral /, eq. (132|) , with the RG scale fi playing the role of the infrared cut-off 
(analogous to that of the system size L in eq. (|34|) ). This yields (in the minimal subtraction 
scheme): 

Z,{t) = l-^^ + 0{t'). (35) 

It is easy to check that this results in a relation between the bare and the renormalized 
coupling constant analogous to eq. 



-fi'- ^ ^ = 1 - -^f.'- . (36) 
t ts ts 2a — 1 
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From eq. (|35D, we get the expression for the /3-function: 



{a-l)t-^ + 0{t'). (37) 



tB 



l + tdt\nZi{t) ' ' 27r 



Both eqs. ( |51D and ( P7| ) show that the couphng constant t increases with the system size L 
(resp. scale /i), which is analogous to the behavior found in the conventional scaling theory 
of localization in < 2 dimensions [^,^ . The RG flow reaches the strong coupling regime 
t ~ 1 at the scale /i ~ t^J^'^ ^\ Remembering the relation of the bare coupling constant ts 
and the index a to the parameters of the original PRBM model: oc 6^""^; a oc 2a — 1, 
we conclude that the length scale 



^ ^ ~ (38) 

plays the role of the localization length for the PRBM model. 

This conclusion is also supported by an inspection of the expression for the IPR, eq. 
(P7|). Evaluating the one-loop perturbative correction in eq. (0), one gets 

= {l + ^(^ - l)^C(-)^^"^} ' (39) 

where ({a) is Riemann's zeta-function; ({c) — l/(o" — 1) for cr close to unity. The cor- 
rection term becomes comparable to the leading (GOE) contribution for the system size 
N ~ t^/^'^~^\ parametrically coinciding with the localization length ^. For larger N the per- 
turbative expression ( P5| ) loses its validity, and the IPR is expected to saturate at a constant 
value {Q ~ ^^^^ for > ^ 

In conclusion of this subsection, let us stress once more that the localized eigenstates in 
the present model are expected to have integrable power-law tails: oc r~^" = r~'^~^ 

at r > ^. 



2q-1 



B. < a < 1 (1/2 < a < 1) 

We start again from considering the perturbative corrections to the diffusion propagator 
(PB|). The integral (^) is now dominated by the region /c ~ g, and is proportional to 
We get, therefore: 



K-\q) = ^-CM (40) 



t 



with a numerical constant C„. 

We see that the correction term is of higher order in |g| as compared to the leading one. 
Thus, it does not lead to a renormalization of the coupling constant t. This is readily seen 
also in the framework of the RG scheme, where the one- loop integral in eq. (|30D gives no 
rise to terms of the form Igl*^. One can check that this feature is not specific to the one-loop 
RG calculation, but holds in higher orders as well. An analogous conclusion was reached for 
the case of a vector model with long-ranged interaction by Brezin et al. [E^ . 
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In our case this implies that the renormahzation constant Z\ is equal to unity and the 
/?— function is trivial: 



m = (1 - . (41) 

This means that the model does not possess a critical point, and, for a < 1, all states 
are delocalized, for any value of the bare coupling constant t. This property should be 
contrasted with the behavior of a d— dimensional conductor described by the conventional 
local non-linear a— model and undergoing an Anderson transition at some critical coupling 



t = tc 21 



Since all states of the model turn out to be delocalized, one might think that their 
statistical properties are the same as for GOE. This is, however, not quite true. In particular, 
calculating the variance of the inverse participation ratio 1^ in the same way as it is done in 
Ref. M we get: 



W-^M!.J^5:nV,/)^j; E^^tfW. (42) 

\ ^1 rr' q=-Kn/N\n=l,2,.. 

At cr > 1 /2 the sum over q is convergent yielding 



Thus, in this regime the fluctuations of the IPR are much stronger than for the GOE where 
S{I) oc 1/A^. Only for ct< 1/2 (a <3/4) the IPR fluctuations acquire the GOE character. 

Considering higher cumulants of the IPR, {{I2)), one finds that the GOE behavior is 
restored at cr < cr^") = 1/n. In this sense, the model is analogous to a d-dimensional 
conductor at d = 2/cr. Therefore, only when a —>■ (correspondingly, a 1/2 in the 
original PRBM formulation) all statistical properties become equivalent to those typical for 
the GOE. 



C. Critical regime: a = 1 (a = 1). 



As we have seen, the case a = 1 separates the regions of localized (a > 1) and extended 
(cr < 1) states. It is then natural to expect some critical properties showing up just at cr = 1. 
Let us again start from considering the generalized diffusion propagator, eq. (|26|) . At a = 1 
the one-loop correction yields 



1 



ITCU 



;K~\q)=A\q\ 



r^-lln(|g|L) 



(44) 



As it was natural to expect for the critical point, the correction to the coupling constant 
is of logarithmic nature. However, eq. (44) differs essentially from that typical for a 2D 
disordered conductor: 



t-' = t-^'-HL/i) , 



(45) 
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where the bare couphng constant ts corresponds to scale /. Comparing the two formulas, we 
see that in eq. (|^) the mean free path / is replaced by the inverse momentum q~^. Therefore, 
the correction to the bare coupling constant is small for low momenta q ~ '^/L, and the 
correlator K{q) is not renormalized. This implies the absence of eigenstate localization, 
in contrast to the 2D diffusive conductor case, where eq. ( ^Sj) results in an exponentially 
large localization length ^ oc expt^^. On a more formal level, the absence of essential 
corrections to the low-g behavior of K{q) is due to the fact that the region k > q does 
not give a logarithmic contribution. This is intimately connected with the absence of t 
renormalization at a < 1. 

To study in more details the structure of critical eigenfunctions, let us consider the set 
of IPR, Iq. The perturbative correction, eq. (pTl) , is evaluated at a = 1 as 



{Q = \^ + QiQ-^){^ In (N/b)] ^^1^ , (46) 



where the microscopic scale b, eq. (|^), enters as the ultraviolet cut-off for the a— model, the 
role usually played by the mean free path /. This formula is valid as long as the correction 

is small: g <^ s^^^ i-^/b) ■ For larger q the perturbation theory breaks down, and one 
has to use the renormalization group approach. This procedure, pioneered by Wegner |17 



and developed by Altshuler, Kravtsov and Lerner p3[ requires the introduction of higher 



vertices of the type Zq J Sti'^ {QkA)dr in the action of the non-linear a— model, and their 
subsequent renormalization. This results in RG equations for the charges Zq which in the 
present case, and in one-loop order, read: 

dZq , t 



q{q-l)—Zq, (47) 



(iln/x~i Stt 
where is the renormalization scale. Integrating eq. (^), we find: 

Note, that this formula is reduced to the perturbative expression, eq. (^61), in the regime 
^ln(iV/6)l"'/'. 



q < 
T 



le behavior described by eq. (H3) is characteristic for a multifractal structure of wave 



functions, when 

(/,) oc iV-^«('?-i) (49) 
with dq being the set of fractal dimensions. We find from eq. 



dq = l-q£. (50) 

This form of the fractal dimensions is similar to that found in two and 2 + e dimensions for 
the usual diffusive conductor [0,ff^,B3]. The one- loop result (EO) holds for q <^ Sn/t. 



The set of fractal dimensions dq (as well as spectral properties at a = 1, see next section) 
is parametrized by the coupling constant t. Strictly speaking, our a— model derivation is 
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justified for t ^ 1 (i. e. b ^ 1). However, the opposite limiting case can be also studied, 
following the ideas of Levitov [0. This corresponds to a (i— dimensional Anderson insulator, 
perturbed by a weak long-range hopping with an amplitude decreasing with distance as 
r^"". The arguments of Levitov suggest that the states delocalize at a < d, carrying 
some fractal properties at cr = rf. Our PRBM model in the limit 6 ^ 1 is just the ID 
version of this problem. This shows that the conclusion about localization (delocalization) 
of eigenstates for a > 1 (resp. a < 1), with a = 1 being a critical point holds irrespective of 
the particular value of the parameter h. Alternatively, the regime of the Anderson insulator 
with weak power-law hopping can be described in the framework of the non-linear a— model, 
eq. (0), by considering the limit t ^ 1. Formally, the non-linear a— model for arbitrary 
t can be derived from a microscopic tight-binding model by allowing n ^ 1 "orbitals" per 
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At any rate, the PRBM model, eqs. (|l|), (^, with arbitrary < 6 < oo, or the a— model, 
eq. (|TB|), with arbitrary coupling constant < t < oo display at cr = 1 a rich critical behavior 
parametrized by the value of h (respectively, t). 



VI. SPECTRAL PROPERTIES. 

Let us consider now the issue of spectral statistics of the PRBM model. As is well known, 
a usual diffusive conductor exhibits the Wigner-Dyson (RMT) level statistics in the limit of 
infinite dimensionless conductance g = 2ttvDL'^~'^ . At finite (7 3> 1, there appear deviations 



2^ , p!9|p5| . In the present section we would like to address the analogous problem in the case 



of the PRBM model. 

The basic quantity characterizing the spectral properties is the two-level correlation 
function 



R{s) 



■{v{E)v{E + uj)) 



(51) 



where s = cu/A, A is the mean level spacing, z/(£') is the density of states at energy and 
(. . .) denotes the ensemble averaging. Following ||19[, we find the leading correction to the 
Wigner-Dyson form E^^{s) of the level correlation function, eq. (|5l|), as 



R{s) 



R 



WD, 



(52) 



where 



c = ^EnV,r') 



1 

iV2 



q=-Kn/N; n=l,2,... 



647r2'^ ' 
t 

const ^. „ „ ^ 



a > 1/2 
a < 1/2 



(53) 



At a < 1/2 the sum divergent at high momenta is cut off at g ~ 7r/6; the procedure leaving 
undetermined a constant of order of unity. 
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The correlation function R{s) is close to its RMT value if o" < 1 (the region of delocalized 
states), or else if cr > 1 and the system size N is much less than the localization length ^, 
eq. p^). Under these conditions, eq. ( ^2] ) holds as long as the correction term is small 
compared to the leading one. This requirement produces the following restriction on the 
frequency s = u/A: 



S < Sc 



a > 1/2 
a < 1/2 



(54) 



At larger frequencies [s > Sc), the form of the level correlation function changes from the 



1/s^ behavior typical for RMT to a completely different one [24|. Extending the calculation 



of Ref. to the present case, we find 



R{s) 



TT 



7rn 
N 



OC 



JYl-l/<7^1/a^l/a-2 ^ 



a > 1/2 ((T ^ 1) 
a < 1/2 



(55) 



At last, let us consider the level statistics in the critical regime a = 1. In this case the 
coefficient of proportionality in the asymptotic expression (|55D vanishes in view of analyticity: 



R{s) 



At 
16^ 



dx 







(56) 



This is similar to what is known to happen in the case of a 2D diffusive conductor p6|,p7 



A more accurate consideration requires taking into account the high-momentum cut-off at 
q ~ h~^. In full analogy with the 2D situation mentioned p6|j28|j27[] , we find then a linear 
term in the level number variance: 



{5N'^{E)) ^ k{N{E)) 

K, = R(s)ds = — . 
J Sir 



(57) 



The presence of the linear term (|57|) (as well as the multifractality of eigenfunctions. Sec. 0) 
makes the case a = 1 similar to the situation on the mobility edge of a disordered conductor 
in d > 2 p8[. Let us finally mention that the value of k, eq. (^7]), is in agreement with the 
formula k = {d — d2)/2d, suggested recently by Kravtsov where d2 is given by eq. 
with q = 2, and d = 1 in the present case. 



VII. NUMERICAL SIMULATIONS 

We have performed numerical simulations of the PRBM model for values of a G [0, 2] 
and b = 1. In Fig. 2 we present typical eigenfunctions for four different regions of a. In 
agreement with the theoretical picture presented above, the eigenstates corresponding to 
a = 0.375 and a = 0.875 are extended, whereas those corresponding to a = 1.25 and 
a = 1.625 are localized. At the same time, one can notice that the states with a = 0.875 
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and a = 1.25 exhibit a quite sparse structure, as opposed to the other two cases. We beheve 
that this can be explained by the proximity of the former two values of a to the critical 
value a = 1.0, where eigenstates should show the multifractal behavior, see Sec. |V[ 

In order to get a more quantitative insight into the properties of the eigenstates, we 
concentrated our attention on the behavior of the mean value of the IPR, (h), and on the 
relative variance, 6 = ((/|) — {h)"^) / {h)"^ ■ At any given a we studied the dependence of 
the quantities {I2) and 5 on the matrix size and approximated these dependencies by 
the power-laws {I2) oc 6 oc 1/A^^ for N ranging from 100 to 2400. On Figs. 3 and 

4 we plotted the values of the exponents u and n obtained in this way, versus the PRBM 
parameter a. The expected theoretical curves following from the results of Sec. [V| are 
presented as well. We see from Fig. 3 that the data show a crossover from the behavior 
typical for extended states {u = 1) to that typical for localized states (z/ = 0), centered 
approximately at the critical point a = 1. We attribute the deviations from the sharp step- 
like theoretical curve to the finite-size effects which are unusually pronounced in the 
PRBM model due to the long-range nature of the off-diagonal coupling. The data for the 
exponent /i (Fig. 4) also show a reasonable agreement with the expected linear crossover, 
/i = 4(1 - a) for 3/4 < a < 1, see eq. (|4|). 

VIII. CONCLUSION. 

In this paper, we have performed a detailed investigation of the REM model with a 
power- law decay of the matrix elements, eq. (^. Physically, one can look at this model 
as describing a particle in a ID disordered system with a power-law hopping term. As 
a theoretical tool, we have used a mapping of the problem onto a supermatrix non-linear 
cr— model, eqs. (pTS]), (|19|). Depending on the value of the power-law exponent a in eq. 
(I) (or, equivalently, cr = 2a — 1 in eq. (p!9[)), three different regimes are found: For 
a > 1 {a > 1) all eigenstates are localized with integrable power-law tails. For a < 1 
(cr < 1) the eigenstates are delocalized, for any value of the bandwidth b of the PRBM 
model, eq. (|) (resp. coupling constant t oc of the cr— model, eq. (|19D). These two 

regimes are separated by the critical value a = 1 (cr = 1), where the structure of eigenstates 
is multifractal, and energy levels show statistics intermediate between Wigner-Dyson and 
Poisson ones. These critical properties are similar to those found on the mobility edge of a 
(i-dimensional disordered conductor. At cr = 1, we find a family of critical points labeled by 
the value of the coupling constant t of the non-linear cr— model, so that the critical behavior 
is parametrized by the value of t. In particular, it determines the multifractal exponents 
dq and the coefficient k, of the linear term in the level number variance, which are given at 
t -C 1 by eqs. ( pOD and (|57|), respectively. 

Turning our attention to the regime of locahzed states, we find that it can be subdivided 
into two domains. In the region a > 3/2 the properties of the model are rather close to 
those of a conventional quasi-lD conductor [Q. On the other hand, for 1 < < 3/2 the 
wave packet spreading on a short time scale is superdiffusive, (|r|) ~ ^1/(2"-!)^ which leads 
to a modification of the Altshuler-Shklovskii "tail" of the spectral correlation function, eq. 
(^), and to an unusual scaling of the localization length ^ oc The regime of extended 

states, a < 1, can be also subdivided into two domains. For a < 1/2, all statistical 
properties are identical to those of the GOE (which corresponds to a = 0). On the other 
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hand, at 1/2 < a < 1 the model is quite similar to a diffusive conductor in d = {a — l/2)~^ 
dimensions. This is reffected, in particular, in the ffuctuations of the IPR (see eq. (|43|) and 
the discussion following it), and in the large-frequency "tail" of the spectral correlator, eq. 

Our conclusion about the existence of a set of critical theories at a = 1, parametrized by 
the coupling constant t, perfectly agrees with earlier results by Levitov [P]. He studied the 
effect of a weak power-law hopping on an Anderson insulator, which corresponds to the limit 
t ^ 1 in the a— model formulation, and arrived at the conclusion of criticality of the model 
at a equal to the spatial dimension d. Let us also note that our results are in accordance 
with the fact of absence of localization effects in the Quantum Fermi Accelerator model [0 . 
As it was pointed out in , the evolution equation for this model takes the form of a finite 
difference equation of the tight binding type with a long range hopping term decaying in a 
power-law fashion, eq. @, with a = 1. Our results show that this case corresponds to the 
critical point with extended eigenstates and intermediate level statistics, in agreement with 
the behavior found in Ref. 0. 

We have presented results of a direct numerical simulation of the PRBM model. Our 
data are in reasonable agreement with the above theoretical picture. However, a more 
detailed numerical investigation of the structure of eigenstates and of spectral statistics is 
certainly desirable. In particular, it would be very interesting to study the critical manifold, 
a = 1, where the multifractal properties of eigenstates and intermediate level statistics are 
predicted by our theory. 
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Figure captions 



Fig. 1 

One-loop diagram contributing to the self-energy part of the density-density correlation 
function. The sohd hne denotes the Q-field propagator, whereas the dashed line represents 
the interaction U{r — r'). 

Fig. 2 

Typical eigenfunctions for the matrix size = 800 and four different values of a: 
&)a = 0.375; h)a = 0.875; c)a = 1.250; d)a = 1.625. 

Fig. 3 

Index u characterizing the dependence of the inverse participation ratio {I2) on the matrix 
size via (72) oc l/N". as a function of a. Points refer to the best-fit values obtained from 
matrix sizes between = 100 and = 1000 (squares) or A^ = 2400 (circles). The dashed 
line is the theoretical prediction for the transition from 1/ = 1, at small a, to = 0, at large 
a. 

Fig. 4 

The same as Fig. 3, but for the index /i, derived from the N dependence of the variance 5 
of the inverse participation ratio: 5 = ((/|) — {12)'^)/ oc The dashed hne corre- 

sponds to the predicted linear crossover from //=latQ;<3/4to// = 0atQ;>l. 
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